Spin-chain system as a tunable simulator of frustrated planar magnetism 
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At low temperatures, weakly coupled spin chains develop a magnetic order that reflects the 
character of gapless spin fluctuations along the chains. Using nuclear magnetic resonance, we identify 
and characterize two ordered states in the gapless region of the antiferromagnetic, Ising-like spin- 
chain system BaCo2V208, both arising from the incommensurate fluctuations along the chains. 
They correspond to the columnar and ferromagnetic ordered states of the frustrated J1-J2 spin 
model on a square lattice, where the spins are encoded in original spin chains. As a result of 
fleld-dependent incommensurate fluctuations and frustrated interchain interaction, Ji can be tuned 
continuously with the magnetic fleld, and its value with respect to a flxed J2 selects the ordered 
state. Spin-chain systems can thus be used as tunable simulators of frustrated planar magnetism. 

PACS numbers: 75.10.Jm, 75.10.Pq, 75.25.-j, 76.60.-k 



The study of important models in condensed-matter 
physics increasingly relies on the concept of quantum sim- 
ulation. This concept is largely associated with the gases 
of cold atoms, whose well-defined and highly tunable 
hamiltonians can be mapped onto many model hamil- 
tonians of interest Cold atoms in optical lattices 

allowed the first simulation of the quantum phase tran- 
sition from a superfiuid to a Mott insulator t^, and, re- 
cently, also the simulation of ordered magnetic states on a 
lattice 0, 13- Scarce difficulties with cold atoms, related 
to the fact that they are gases, are absent in quantum 
antiferromagnets, which also feature well-defined hamil- 
tonians. Among them, those that contain dimers of quan- 
tum spins can host bosonic excitations. As the density 
of bosons can simply be driven by the applied magnetic 
field, these systems have proven as convenient simula- 
tors of Bose-Einstein condensates [Hl-Q and Tomonaga- 
Luttinger liquids (TLLs) Q. 

In exploring the field-temperature phase diagram 
of the antiferromagnetic, Ising-like spin-chain system 
BaCo2V2 08 tS-lSj, we discovered a new low-temperature 
ordered state above Bp = 8.6 T. The state borders on the 
recently identified, novel incommensurate (IC) ordered 
state between Be = 3.9 T [3j 15 1 and Bp, shares its 
quantum origin, but exhibits a different order. We show 
that if we map the configuration of spins in each chain 
onto a classical spin, the two ordered states correspond 
to the columnar and ferromagnetic ordered states of the 
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from ferromagnetic to antiferromagnetic, and we extract 
this field dependence. The value of Ji with respect to 
a fixed exchange coupling J2 between the next-nearest- 
neighboring (NNN) chains (i.e., classical spins) then se- 
lects the most favorable ordered state [20|. This is a 
mechanism for the transition between the two ordered 
states at Bp. From the practical point of view, the ex- 
change coupling scheme encountered in BaCo2V208 al- 
lows a simple quantitative control of the exchange interac- 
tion between the effective classical spins, extending the 
simulation possibilities of quantum antiferromagnets to 
the domain of planar magnetism. This is all the more 
valuable as a direct control of exchange interactions in 
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frustrated J1-J2 spin model on a square lattice 
Remarkably, the effective exchange coupling Ji between 
the nearest-neighboring (NN) chains (i.e., classical spins) 
can be continuously tuned by the applied magnetic field 



known materials, by changing the applied pressure 
or the chemical composition [22], is only qualitative. 

We start by explaining the principle of the tuning 
mechanism. It is based on two ingredients: (i) chains 
hosting a type of spin fluctuations that can be manipu- 
lated by the magnetic field, and (ii) a frustrated interac- 
tion between pairs of chains. The tuning mechanism ap- 
plies to chains formed by S" = 1/2 spins that are coupled 
by an anisotro pic, Ising-like antiferromagnetic exchange 
interaction '23'-'2d| . When the magnetic field points along 
the exchange-anisotropy axis z, a single chain exhibits 
gapless spin fluctuations in a certain field range, thus re- 
alizing a TLL [23|. In this range, i.e., between the critical 
fields Be and JSg, a uniform magnetization per spin 
induced by the field B grows monotonously with B from 
to the saturated value of 1/2. In the lower part of this 
range, gapless fluctuations are longitudinal (i.e., they in- 
volve spin components along the fleld) and appear at the 
field-dependent IC wavevector 2kF{B) = tt[1 — 2mz{B)] 
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FIG. 1: (color online) (a),(b) A frustrated pair of antifer- 
romagnetic Jif couplings in the zigzag scheme sums up to 
the effective coupling Ji{B) given by Eq. ([T]), if the chains 
are hosting ordered SDWs with the wavevector 2kF(B). (c) 
Ji{B) based on the numerically calculated m^i^B) for ex- 
change anisotropy A = 2.17 in the range between the critical 
fields Be and Bs- SDWs are mapped onto the classical spins 
(arrows) that assume a ferromagnetic (antiferromagnetic) ar- 
rangement for a negative (positive) Ji{B). 



along the chain [2^ . Mapping the spin chain onto a one- 
dimensional (ID) system of spinless fermions using the 
Jordan- Wigner transformation 27], this equation sim- 
ply relates the filling of the fermion band 1/2 — with 
the wavevector 2kp connecting both points of the Fermi 
surface. Gapless fluctuations prevent an isolated chain 
from ordering even at zero temperature. However, in 
weakly coupled chains, a 3D magnetic order develops at 
low enough temperatures, whose type reflects the char- 
acter of fluctuations: spins develop an ordered compo- 
nent along the field (on top of the uniform niz) in the 
form of spin-density waves (SDWs) with an IC wavevec- 
tor 2kp{B) along the chains [28]. The main point here is 
that the wavevector of SDWs is controlled by the field, 
varying from tt to between the critical fields. This is 
our first ingredient. 

To make the energy of the exchange-coupled SDWs 
depend on their wavevector, and hence on the magnetic 
field, pairs of SDWs should interact in a frustrated man- 
ner. This is our second ingredient. Namely, in the case 
of a non-frustrated interaction based on the couplings Ji 
[Fig.[T][b)], the energy of the coupled pair of SDWs is pro- 
portional to ^ Ji cos((/?i — 1P2), where ipi and (p2 are their 
phases. There is no wavevector dependence. In contrast, 
for the simplest case of a frustrated interaction based on a 
zigzag scheme with the couplings Ji/ [Fig.[lja)], the cou- 
pling energy is proportional to Ji/ cos(2fcF) cos(<y9i — V2). 
The wavevector-dependent cos{2kp) factor simply re- 
fiects a one-site offset between the coupled sites in the 
respective chains. A comparison of the cases displayed 
in Figs. [TJa) and (b) shows that each frustrated pair of 
Jif couplings sums up to the field-dependent effective 
coupling 



Furthermore, due to the cos((/3i — (P2) factor the energy 
of the coupled SDWs adopts the standard form of the 
coupling energy between the classical planar spins with 
orientations (pi and (^2 ■ The phases of SDWs thus encode 
the orientations of the effective classical spins interacting 
with the coupling Ji{B) given by Eq. ([T]). This coupling 
can be tuned continuously by the field B in the range 
between the critical fields, where it even changes sign 
[Fig. Hie)], which is an unexpected feature. Although 
the individual interchain couplings Ji / are antiferromag- 
netic, Ji{B) can be tuned from ferromagnetic to antifer- 
romagnetic. Accordingly, the effective classical spins flip 
from a parallel to an antiparallel arrangement [Fig.[TJc)], 
implying that the original SDWs switch from being in 
phase to being in opposite phase. 

In the following we demonstrate that such a tuning 
mechanism is realized in BaCo2V2 08, which contains 



parallel chains of magnetic Co + ions. 



each carrying 
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an effective anisotropic 5' = 1/2 spin 
chains run along the c axis and form a square lattice in 
the tetragonal a-h plane [Figs. [2fa)-(f)]. The exchange- 
anisotropy axis coincides with the chain direction, which 
is also chosen as the field direction. Interaction between 
the NNN chains (neighbors along the bisectors of a and 
6) is not frustrated and is given by the field-independent 
coupling J2 [Figs, [lib ) and (e)]. A coupling scheme be- 
tween the NN chains (neighbors along a or h) is frus- 
trated, but more complicated than the simple zigzag 
scheme. It contains the couplings of four different sizes 
[expressed by the weights joi ji and j2 relative to J2 
in Fig. [2l^d)], which combine into eight frustrated pairs 
within four chain units [Fig. [Hd)]. As each pair con- 
tributes a term similar to Eq. ([1]), together they amount 
to the effective coupling 

Ji{B) = 4J2 [io+(j-i+Ji) cos[2fcF(S)]+j2 cos[4fcF(S)] 

(2) 

The non-monotonous shape of Ji{B) displayed in 
Fig. Eljg) is generic, while the precise values depend on 
jo, j-i +ji and j2. 

According to Fig. [^f), in the picture with SDWs 
mapped onto the classical spins, a coupled chain sys- 
tem in BaCo2V208 maps onto the frustrated J1-J2 spin 
model on a square lattice. NN spins in this model interact 
with the tunable coupling Ji (B) and NNN spins with the 
field-independent coupling J2. These two couplings com- 
pete to determine the most favorable spin arrangement. 
The effective classical spins are free to rotate in a plane 
and, for the arrangement described by the wavevector 
k = {kx,ky), the energy of the model reads [2^ 
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cos k„ 



J9 cos k^ cos k,. 



(3) 



Ji{B) = 2Ji/cos[2fcF(B)]. 



(1) 



Depending on the relative values of Ji{B) and J2, this en- 
ergy is minimized in one of the three ground states: ferro- 
magnetic with k = for Ji{B) < — 2J2 [Fig.[21[i)], colum- 
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FIG. 2: (color online) (a)-(c) Crystal structure of BaCo2 V20g contains edge-sharing CoOe octahedra forming the chains along 
c. The dominant exchange interaction between Co^^ spins runs along the chains. Much smaller exchange interaction between 
spins in neighboring chains is mediated by the O-O bridges. Corresponding schemes of interchain exchange couplings are 
shown in (d)-(f). (d) Couplings of four different sizes [highlighted, also shown in (a)] between the NN chains combine into eight 
frustrated pairs, which sum up to the tunable effective coupling Ji (B) . (e) Non- frustrated interaction between the NNN chains 
is based on couplings J2. (f) Chains form a frustrated square lattice in the a-b plane, where Ji{B) competes with J2. (g) Ji{B) 
given by Eq. ([2| calculated for jo = 1.0, j-i + ji = 3.4, j2 = 2.3 [all defined in (d)] and exchange anisotropy A = 2.17 between 
the critical fields Be ~ 3.9 T and Bs — 22.7 T. The value of Ji{B) with respect to J2 selects the ordered state. (h),(i) Sketches 
of the columnar and ferromagnetic ordered states with SDWs mapped onto the classical spins (arrows), (j) Field-temperature 
phase diagram obtained from the NMR data shown in Figs. [4{a) and (b) for a field parallel to c. It agrees with the specific 
heat data available up to ~_Bp — 8.6 T 14]. Lines are calculated phase boundaries based on Ji{B) [plotted in (g)] between the 
TLL state and the ordered states. 



nar with k = (7r,0) or (0,7r) for -2J2 < Ji{B) < 2J2 
[Fig. [2jh)], or checkerboard with k — (7r,7r) for Ji{B) > 
2J2 (not shown) In the original picture, k describes 
the arrangement of the phase shifts of SDWs in the a-b 
plane. Finally, as we are free to tune Ji{B) in accordance 
with Fig. [DJg), covering all three ground states, we are 
in possession of a simulator of frustrated magnetism on 
a square lattice. 

We check this picture experimentally by extracting 
the spin arrangement in various states from the corre- 
sponding ^^V nuclear magnetic resonance (NMR) spec- 
tra [2^. At temperatures above the ordered states, 
throughout the TLL state, the spectrum exhibits a sin- 
gle peak [Fig.[3][a)] generated by a single crystallographic 
V site located in between the NN chains [Figs. Ufa) and 
(c)]. This symmetry breaks at the transition to one of the 
ordered states, and the spectrum develops features. The 
appearance of the sinusoidal SDW breaks the transla- 
tional symmetry along the chains and leads to the charac- 
teristic U-shaped spectrum [sOl. Between Be = 3.9 T and 
Bp = 8.6 T, the spectrum develops two U-shaped com- 
ponents [Fig. |31Ib)], a fingerprint of the columnar SDW 



state. The components are generated by two sets of V 
sites arising from the broken symmetry between the a 
and h directions. Above Bp and up to at least 15 T, the 
highest field in our experiment, the spectrum consists of a 
distribution of U-shaped components, which are revealed 
during the transverse (i.e., T2) relaxation [Figs.[3][c) and 
(d)]. The distribution is generated by a single set of V 
sites, implying that the symmetry between a and h direc- 
tions is restored. This suggests that either the ferromag- 
netic or the checkerboard SDW state is realized above 
Bp. Small linewidth of the spectrum is compatible with 
the ferromagnetic SDW state [2§|. The spectrum is per- 
fectly reproduced by a distribution that corresponds to 
the 2D sinusoidal modulation of the SDW amplitude in 
the a-b plane [Fig. |31[c)]. 

The field dependence of Ji is deduced from the mea- 
sured shape of the phase boundary between the ordered 
SDW states and the TLL state [Fig. WL})]- Transition 
temperatures in different fields are obtained from the 
temperature dependence of the '^^V longitudinal and 
transverse relaxation rates, T^^ and T^^ [29| . In par- 
ticular, the divergent behavior of T^^(T) and T2^{T) on 
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FIG. 3: (color online) Representative ^^V NMR spectra of 
various states (for selected field and temperature values) are 
recorded in a field parallel to c after different durations of T2 
relaxation (given), (a) In the TLL state the spectrum exhibits 
a single peak, (b) In the columnar SDW state, the unrelaxed 
spectrum is reproduced by two characteristic U-shaped com- 
ponents of the same area. They are generated by two sets 
of V sites, located in between the chains hosting SDWs that 
are either in phase (dashed line) or in opposite phase (dotted 
line). The T2 relaxation preserves the shape of the spectrum, 
(c) In the ferromagnetic SDW state, the unrelaxed spectrum 
is reproduced by a distribution of U-shaped components gen- 
erated by a single set of V sites (dashed line). A distribution 
based on the 2D sinusoidal modulation of the SDW ampli- 
tude in the a-h plane results in a spectrum with a large cen- 
tral weight and pronounced wings (indicated by arrows). A 
longer T-i, relaxation suppresses the narrow U-shaped compo- 
nents and reveals the broad ones, (d) Broad components from 
(c) scaled up 10-times. 



Whereas the field- independent J'{B) below Bp implies 
a negative slo pe of the boundary (already known be- 
fore our work (l4i [28|). the slope inversion above Bp is 
due to J'{B) starting to grow steeply enough with field 
due to an increasing ~Ji{B) [Fig. ^g)]. As the mean- 
field expression for the transition temperature assumes 
a uniform SDW amplitude in the a-b plane, a perfect 
agreement with the experimental points suggests that the 
observed 2D amplitude modulation in the ferromagnetic 
SDW state is nearly uniform, i.e., is long wavelength. We 
also note that the extracted Ji{B) in Fig. ^g) predicts 
the existence of the checkerboard SDW state above 20 T, 
which remains to be studied in high-field experiments. 

Finally, we verify that the behavior of individual spin 
chains in the covered field range is indeed dominated by 
the longitudinal IC spin fluctuations, one of the ingre- 
dients for the concept of tunable interaction. Gapless 
fluctuations are directly probed by Tf ^ At 2.4 K, 

in the TLL state, the Tj~^ data precisely follow the de- 
creasing trend of the TLL prediction for the longitudi- 
nal IC fluctuations [Fig. IH^c)]. There are no traces of 
the transverse staggered fluctuations, which are expected 
to become important at higher fields 28|. The neces- 
sary condition for the dominance of the longitudinal IC 
fluctuations is the anisotropy of the involved S = 1/2 
spins [2^. An effective description of this type indeed 
applies to Co^"*" spins in trigonally distorted octahedral 
environments the case encountered in BaCo2V208. 
Related Co-based spin-chain systems show promise for 
use as simulators of frustrated magnetism in other planar 
geometries, such as triangular or kagome, hence allowing 
new insights into spin-liquid physics [33 1. 



increasing temperature T is a signature of the second- 
order phase transition [Figs. SJa) and (b)]. The phase 
boundary changes the slope at exactly Bp, on the bor- 
der between the two SDW states. In the columnar SDW 
state, Eq. ^ yields the coupling energy e oc —J2, while 
in the ferromagnetic SDW state it is £ oc Ji (B) + J2 . The 
average coupling of the chain to its neighbors in the two 
ordered SDW states is thus 



J'{B) 



J2 if B < Bp 

-MB) -J2 \iB>Bp 



(4) 



which is continuous at Bp because Ji{Bp) — — 2J2. For 
weakly interacting chains, the interchain coupling J' can 
be treated in the mean-fleld approximation [Ssl 31 1 . 



This basically implies that the transition temperature 
scales as a field-dependent power of J'. Using J'{B) 
given by Eq. (j4|), the measured phase boundary is per- 
fectly reproduced with the choice of Ji{B) plotted in 
Fig.lHg) and the interchain coupling of J2//CB = 41 mK 
(ks is the Boltzmann constant), which is indeed much 
smaller than the intrachain coupling of 65 K 
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In summary, we identified a mechanism for the tun- 
able exchange interaction between the NN spin chains 
in BaCo2V208- This mechanism is responsible for the 
transition to a new IC ordered state that we discovered, 
and allows to regard BaCo2V2 08 as a tunable simulator 
of frustrated magnetism on a square lattice. The novel 
phenomenon of tunable interaction emerges in the pairs 
of ID objects dominated by tunable quantum fluctua- 
tions in the presence of frustration. It adds to the list of 
other emergent phenomena arising from particular types 
of spin fluctuations in frustrated spin systems, such as 
emergent gauge fields in classical spin ices or fractional 
particle excitations in quantum spin liquids [33| . 
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FIG. 4: (color online) (a),(b) ^^V longitudinal and transverse relaxation rates, T^^ and T^^, as a function of temperature T 
in different magnetic fields parallel to c. Arrows mark the extracted temperatures of transition from the ordered SDW state 
to the TLL state. Solid lines are guides for the eye. (c) Interpolated T^^ from (a) as a function of the field B at 2.4 K (in the 
TLL state) is compared to the TLL prediction for two types of spin fluctuations. 
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SUPPLEMENTAL MATERIAL 

Experimental details. We grew a high-quality sin- 
gle crystal of BaCo2V208, shaped as a 5 cm long rod 
of 3 mm diameter, using the floating zone technique in a 
mirror furnace under air. Prom the rod, a sample with di- 
mensions 2x0.8x0.2 mm"'* was cut, so that the c axis was 
perpendicular to the biggest face. Nuclear magnetic reso- 
nance (NMR) experiments were performed on ^^V nuclei 
with spin / = 7/2 in a magnet that enabled field sweeps 
up to 15 T. A copper NMR coil with the sample was in- 
stalled into the mixing chamber of a dilution refrigerator, 
so that the c axis of the sample was parallel to the field. 
^^V spectra were measured by sweeping the frequency 
at a fixed magnetic field and summing the Fourier trans- 
forms of the echo. To be able to cover wide frequency and 
field ranges, we used the "top tuning" scheme, where a 
variable tuning delay line and a matching capacitor were 
mounted outside the NMR probe. The obtained NA^IR 
spectra were corrected for the demagnetizing effect by 
means of the known static magnetization curve. The field 
was calibrated by ^''Cu NMR of metallic copper in the 
NMR coil, and the temperature was measured using a 



calibrated Ru02 resistor. The longitudinal (i.e., Ti) re- 
laxation curves were obtained from the integrated spin- 
echo intensity after the variable time delay following an 
inversion pulse. Due to the small quadrupole splitting, 
amounting to 30 kHz, all nuclear transitions could be ir- 
radiated at once, so that the relaxation curves could be 
reproduced by a simple exponential function. In the or- 
dered phase, we needed to add the stretching exponent. 
The transverse (i.e., T^) relaxation curves were obtained 
from the integrated spin-echo intensity with variable time 
delay between the two pulses. Characteristic quadrupole 
oscillations prevented us from fitting the curves with the 
standard decay functions, so that the T2 values were 
estimated as the delays at which the intensity reached 
1/e = 0.37 of its maximum value. The details of mag- 
netic ordering were extracted from the ^-'^V NMR spectra 
using the tensors of the hyperfine and dipolar couplings 
between the ^^V nucleus and the Co^+ spins. The hy- 
perfine tensors were determined on the basis of our high- 
temperature study of orientation-dependent ^^V NMR 
spectrum, while the dipolar tensors were calculated nu- 
merically from the known crystal structure. 



